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Abstract 

We investigate the B meson light-cone distribution amplitudes in the heavy- 
quark limit which are relevant for the QCD factorization approach for 
the exclusive B meson decays. We derive exact relations between two- 
and three-particle distribution amplitudes from the QCD equations of mo- 
tion and heavy-quark symmetry constraint. As solution of these relations, 
we give representations for the quark-antiquark distribution amplitudes in 
terms of independent dynamical degrees of freedom. In particular, we find 
that the Wandzura-Wilczek-type contributions are determined uniquely in 
analytic form in terms of A, a fundamental mass parameter of heavy-quark 
effective theory, and that both leading- and higher-twist distribution am- 
plitudes receive the contributions of multi-particle states with additional 
gluons. 



*JSPS Research Fellow. 



Recently systematic methods based on the QCD factorization have been developed 
for the exclusive B meson decays into light mesons |Ij, @, |J §] (for other approaches 
see, e.g. ||). Essential ingredients in this approach are the light-cone distribution am- 
plitudes for the participating mesons, which constitute nonperturbative long-distance 
contribution to the factorized amplitudes 0. As for the light mesons (ir, K, r], p, u, 
K*, 0) appearing in the final state, systematic model- independent study of the light- 
cone distributions exists for both leading and higher twists || || [K], [TTJ] . On the other 
hand, unfortunately, the light-cone distribution amplitudes for the B meson are not 
well-known at present and provide a major source of uncertainty in the calculations of 
the decay rates. In this Letter, we demonstrate that heavy-quark symmetry and con- 
straints from the equations of motion determine a unique analytic solution for the B 
meson light-cone distribution amplitudes within the two-particle Fock states. We also 
derive the exact integral representations for the effects of higher Fock states with addi- 
tional gluons. A complete set of the B meson distribution amplitudes are constructed 
in terms of independent dynamical degrees of freedom, which satisfies all relevant QCD 
constraints. 

In the heavy-quark limit, the B meson matrix elements obey the heavy-quark 
symmetry, and is conveniently described by the heavy-quark effective theory (HQET) 
13| . Following Refs. |14], H, U, we introduce the quark-antiquark light-cone distri- 



bution amplitudes <fi±{t) of the B meson in terms of vacuum-to-meson matrix element 
of nonlocal light-cone operators in the HQET: 



(01^(^(0)1^)) = -^^ 



r 1+ 111 m ,<f>+(t)-4>-(t) 

75r^— < (f>+(t)- i— 



(i) 



2t 

where z 2 = 0, v 2 = 1, t = v ■ z, and = Mv^ is the 4-momentum of the B meson with 
mass M. h v (x) denotes the effective 6-quark field, b(x) ~ exp(— im^v ■ x)h v (x), and is 
subject to the on-shell constraint, tfh v = h v |T^, [Tj|. T is a generic Dirac matrix and, 
here and in the following, the path-ordered gauge factors are implied in between the 
constituent fields. Jb is the decay constant defined as usually as 

(0|?(0)y*75^(0)|S(p)) = i/ B M/ , (2) 

so that <p±(t = 0) = 1. The behavior of the RHS of eq.(|l]) for a fast-moving meson, 
t = v ■ z — >• oo, shows that 0+ is the leading-twist distribution amplitude, whereas 0_ 
has subleading twist |14 |. 

It is well-known that the equations of motion impose a set of relations between 
distribution amplitudes f| [| [l(| |TT1 . To derive these relations, the most convenient 
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method is to start with the exact operator identities between the nonlocal operators: 
q{x)^Th v (0) = q{x) p Th v {0) + i I duu q{x)gG^(ux)x u ^Th v (ti) , (3) 



o 



dx^ 

d _ r 1 

— — q(x)Th v (0) = —q(x)TDnh v (Q)+i I du(u — 1) q(x)gGn U (ux)x L Th v (0) 
(Jxt 1 Jo 

+ ^{g(x)r^(0)} , (4) 

where x^ is not restricted on the light-cone. _D M = <9 M — igA^, D^—dfj, +igA fl are the 
covariant derivatives, G^ u = {if ' g)[D^D v ] is the gluon field strength tensor, and 



d„{q(x)rh v (0)} = q(x + y)Th v (y) 



(5) 

y^O 



stands for the derivative over the total translation. These identities simply describe the 
response of the nonlocal operators to the change of the interquark separation and/or 
total translation. 

By taking the vacuum-to-meson matrix element of these relations, the operators 
involving the equations of motion vanish. Going over to the light-cone limit x M — ► z^, 
the remaining terms can be expressed by the appropriate light-cone distribution ampli- 
tudes (f> + , </>_, etc. The terms given by an integral of quark- ant iquark-gluon operator 
are expressed by the three-particle distribution amplitudes corresponding to the higher- 
Fock components of the meson wave function. Through the Lorentz decomposition of 
the three-particle light-cone matrix element, we define the four functions tyy (t , u), 
a (t , u), X A {t , u) and Ya (t , u) as the independent three-particle distributions: 

(0\q(z) gG, u (uz) z»Th v (0)\B(p)) 
= l -f B MTv 75 r^|(^/-t 7/1 ) (* A {t,u)-* v {t,u)) (6) 

-ia^z v ^ v {t,u)- ZpX A (t,u) +^ jY A (t,u)\ . 

This is the most general parameterization compatible with Lorentz invariance and the 
heavy-quark limit. 

By using q p= for the light quark and substituting eqs.(Q) and (|B|), the first iden- 
tity (f|) yields the two differential equations connecting the two-particle distributions 
<p + and 0_ with the three-particle distribution amplitudes: 

4>'_(t) | (0+(*) — ^-(*)) =2*^ duu(* A (t,u)-* v (t,u)) , (7) 
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4>'+(t) - ~^-{t)-\{M)-^{t))+At d -^ 

= 2t f duu (V A (t, u) + 2 v (t, u) + Xi(f , u)) , (8) 



(9) 



where <p±(t) = d<f)±(t)/dt, and we introduced a shorthand notation 

90+ (t) _ d(f) + (t,x 2 ) 



dz 2 <9x 2 



Here, via 0+(t) — > 4>+(t,x 2 ), we extend the definitions in eq. (HD to the case z — > x 
(a; 2 7^ 0), since the derivative in the LHS of eq.@ has to be taken before going to the 
light-cone limit. 

Next, we proceed to the second identity To use the HQET equation of motion 
v ■ Dh v = |13| , we contract the both sides of eq.(|]) with v^. After similar 
manipulations as above, we finally obtain another set of two differential equations: 

1 (1 ^ 1 U<\ , -XI , nM+W 



4>'+(t) -(4> + (t)-4>4tj)+iA4>+(t) + 2t- 



2t > ^ v ') dz 2 

tf du(u-l)(V A {t,u)+X A {t,u)) , (10) 







ft + (t) i'_(t) + (iA--}(i + (t)-4>-(tj)+2t 



dz 2 dz 2 



2t [ du {u - 1) $ A {t, u) + Y A {t, u)) , (li; 
J o 



where 



iv d^\qTK\B{p)) . . 

A = M — rrih = — r~, ; — , - , N , (12) 

(0\qTh v \B(p)) 1 ; 

is the usual "effective mass" of meson states in the HQET |15|, 16 , 12 . 

Eqs.(|7D- (|TTD are exact in QCD in the heavy-quark limit, and are the new results. 

In the approximation that the three-particle amplitudes in the RHS are set to zero 

( "Wandzura-Wilczek approximation"), the system of differential equations (0), is 

equivalent to that found in Ref. 0. By combining eqs.(|) and (PUD, we can eliminate 

the term dcf) + {t)/dz 2 as 

4>'+{t) + #_(*) +2iA$ + (t) = -2t£du (v A (t,u) +X A (t,u) +2utf v (t,ii)) . (13) 

Now, from a system of equations (|?p and ([T3"D, we can obtain the solution for + and 
</>_. For this purpose, it is convenient to work with the momentum-space distribution 
amplitudes <j)±{uS) defined by 

0±(t) = f du e~ iwt (j) ± {u) . (14) 



3 



Here uv + has the meaning of the light-cone projection k + of the light-antiquark mo- 
mentum in the B meson. The singularities in the complex t-plane are such that <j)±(u)) 
vanish for u < 0. Similarly for the three-particle distribution amplitudes, we define 

F(t, u)= J dudi e-^^Fiu, , (F = {V v , X a , Y a }) . (15) 

Here uv + and £v + denote the light-cone projection of the momentum carried by the 
light antiquark and the gluon, respectively, and F(u,£) vanishes unless u > and 
£ > 0. By going over to the momentum space, we obtain from eqs.(^) and (p~3|), 

+ + M = /H, (16) 

duj 

(u - 2A) <j>+(u) + u<p_(uj) = J{uj) , (17) 
where I(u) and J{oS) denote the "source" terms due to three-particle amplitudes as 

= 2 ^J??£M 1 * aM -* vM] - (18) 
jh = -2-f f"d P r ^[* A (p,()+x A (p,d] 

dui Jo Jlu-p £ 

d j- |9^e (19) 

We solve eqs.([T6|) and ( |17D with boundary conditions 0±(o>) = for c<j < or uj — >■ oo. 
Obviously, the solution can be decomposed into two pieces as 

± M = 0i w) M + #H , (20) 

where (j)^ W \u) are the solution of eqs.([l6|) and fllTD with the source terms set to zero, 
I(uj) = J(uj) = 0, i.e., in the Wandzura-Wilczek approximation. <fr±\u)) denote the 
pieces induced by the source terms. 

First, let us discuss (jy± W . Eq.([L6]) alone, with I(u) = 0, is equivalent to a usual 
Wandzura-Wilczek type relation derived in Ref. [0: 

fdp^W. (21) 

Now combining eqs. (fi6|) and (17) , we are able to obtain the analytic solution explicitly 

as 

rffV) = ^»(2A--), (22) 

4>r vi n = ^ fl(2 A- w) , (23) 
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for u > 0. These are normalized as f£° du(jy± (ct>) = 1. They vanish for uo > 2A, 
and we note that 2A is actually the kinematical upper bound of u allowed for the 
two-particle Fock states of the B meson in the heavy quark limit. 
The solution for can be obtained straightforwardly and reads: 

= (24) 

&\») = 2 A^H-) + m , (25) 



2A v ' uo 



where u > and 



with 



ifW=/W + (^-^)jW. (27) 

These functions obey Jq°° duj<$± (uj) = 0, so that the total amplitudes are normalized 
as Jq°° dui(j>±{u}) = 0±(O) = 1. The solution Q2Tj) with eqs.fl2"2])-(P7D for the distribution 
amplitudes is exact and presents the principal result of our paper. Eq.(^) might ex- 
hibit discontinuity at u = 2A, but this does not constitute a problem because physical 
observables are given by convolution integrals of distribution amplitudes with smooth 
coefficient functions, i.e., the sharp behavior is averaged over. Namely, the distribu- 
tion amplitudes should generally be understood as distributions (in the mathematical 
sense) . 

In terms of the Mellin moments (n = 0, 1, 2, ■ • •), 

(u n )± = I duu n M^) = / du u n (i) { ± W \uj)+ du uj n (f) { l\uj) 

_ {ujn) (WW) + {ujn){ g) ^ (2g) 







our solution reads: 



W) = -^(2A)" , = ( — y, -^ (2A)" , (29) 

n + 2 (n + l)(n + 2) 



o n— 1 




5 



+ ( n + 2-i)[X A ]^ i + (n + 3-i)^-[* v }]- i ) , (30) 



= 4r^^^£r i^Nr-ini' ( 3i ) 

n + 1 ri + 1 ~[ \ J/J + 1 V 7 

where (^j = — and we introduced the double moments of the three-particle 

distributions as 

poo poo 

[F]) = J o du j o d£ J-^^F(u, £) , (F = {^ v , ¥ A , X A }) . (32) 



For a few low moments, eqs.(Pq)-(|32|) give 



M+ = ~A, (">- = ~A, (33) 

<^ 2 >+ = 2A 2 + ^A| + ^, ^)- = \a 2 + \\ 2 h . (34) 
Here we have used 

M; = ^A|, [*v]i = iA^, [X A ]1 = 0, (35) 

where A# and Xh parameterize the two independent reduced matrix elements of local 
quark- ant iquark-gluon operators of dimension 5, and are related to the chromoelectric 



and chromomagnetic fields in the B meson rest frame |14 



(0\qgE-cx l5 h v \B(p=0)) = f B M\% , (36) 
(0 1 qgH- a l5 h v \B(p = 0)) = if B M\ 2 H , (37) 

with E i = G 0i , H i = -\e^ k & k , and ct = 7 °7. The results (||) and (0) exactly 
coincide with the relations obtained by Grozin and Neubert ||14j| , who have derived 
their relations by analyzing matrix elements of some local operators. Our results fl28|)- 
([32|) from nonlocal operators give generalization of theirs to n > 3. 

An interesting feature revealed by our results is that the leading-twist distribution 
amplitude <p + as well as the higher-twist 0_ contains the three-particle contributions. 
This is in contrast with the case of the light mesons, where the leading-twist amplitudes 
correspond to the "valence" Fock component of the wave function, while the higher- 
twist amplitudes involve contributions of multi-particle states. 

In the higher-twist distribution amplitudes of light mesons HI, the contri- 

butions of multi-particle states with additional gluons have been generally important 
and broadened the distributions, but they have typically produced corrections less 
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than ~ 20% to the main term given by the Wandzura-Wilczek contributions. We 
note that, in the present case, there exists a rough estimate A^/A 2 = 0.36 ± 0.20, 
A^/ A 2 = 0.60 ± 0.23 by QCD sum rules |TJ (see eq.(0)), but any estimate of the 
higher moments is not known. It is obvious that further investigations are required 
to clarify the effects of multi-particle states. In the applications to the physical am- 
plitudes, the evolution effects including the three-body operators also enter the game. 
All these further developments for going beyond the Wandzura-Wilczek approximation 
can be exploited systematically starting from the exact results in this paper, as it has 
been done for light mesons. 

Inspired by the QCD sum rule estimates, Grozin and Neubert 0] have proposed 

GN I . ,\ _ / w \ ( u \ jJGNf, ,\ _ f 1 



model distribution amplitudes (jy^_ (u>) = ypzj exp \— ^J, 0_ (u>) = y^j exp ^— ^ 
where ujq = 2 A/3. The shape of their model distributions is rather different from that of 
the Wandzura-Wilczek contributions ( p2]) and (|23|), except the behavior <p9 N (uj) ~ u, 
(f)_{uj) ~ const, as uo — ■> 0. Actually, such behavior when the light-antiquark becomes 
"soft" is suggested by the corresponding behavior of the light-mesons || |9|, [TO], [XT 



But we note that the gluon correction (|24| ) to 4>+{<~o) would modify such behavior if 

J(0) = -2/ °°(^/0 (^(0,0 + ^(0,0) ^ 0. Q 

To summarize, the solution in this paper provides the powerful framework for build- 
ing up the B meson light-cone distribution amplitudes and their phenomenological ap- 
plications. Our results represent the quark-antiquark distribution amplitudes in terms 
of independent dynamical degrees of freedom, and satisfy the constraints from the 
equations of motion exactly. The other essential constraints are imposed by heavy- 
quark symmetry, which lead to the reduced set of the distribution amplitudes and 
thus allow us to determine the Wandzura-Wilczek contributions explicitly in analytic 
form. The Wandzura-Wilczek contributions give the effects corresponding to the va- 
lence distributions, while the corrections due to the additional gluons are also obtained 
as exact integral representations involving the three-particle distribution amplitudes. 
A detailed study of the gluon corrections requires systematic treatment of the three- 
particle distributions, and will be presented elsewhere. 

We have treated the light-cone distribution amplitudes for the pseudoscalar B me- 
son. In the heavy-quark limit, the distribution amplitudes of the vector meson B* is 
related with those of the B meson thanks to heavy-quark spin symmetry |TJ| , so that 
the solution in this paper also determines a complete set of the B* meson distribution 
*J.g)t 



(0) of eq.(E3) is finite even if J(0) ^ 0. 
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amplitudes. 

Among the four differential equations (pj)-(|TT|) derived in this paper, we have utilized 
only a system of two equations to obtain the solution (pUP- Now, with our explicit 
solution, the other two equations can be used to determine d(f>±{t) / dz 2 (see eq.©), 
i.e., the transverse momentum dependence of the distributions, which are necessary for 
computing the power corrections to the exclusive amplitudes. 

The authors would like to thank T. Onogi for fruitful discussions. The work of J.K. 
was supported in part by the Monbu-kagaku-sho Grant-in- Aid for Scientific Research 
No.C-13640289. The work of C-F.Q. was supported by the Grant-in-Aid of JSPS 
committee. 
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Erratum to: "B meson light-cone distribution amplitudes 

in the heavy-quark limit" 
[Phys. Lett. B 523 (2001) 111] 

Hiroyuki Kawamura a , Jiro Kodaira b , Cong-Feng Qiao b , Kazuhiro Tanaka c 

a DESY-Zeuthen, Germany 
b Department of Physics, Hiroshima University, Japan 
c Department of Physics, Juntendo University, Japan 



There was an error in Eq. (6) of our Letter. The correct equation is 



{0\q(z)gG^(uz) z v T h v (0)\B(p)) 



7 5 r^i|(^ t-t llx ) (v A (t,u)-y v (t, 



u 



(6) 



i a^z v iffv(t, u) - Zf, X A (t, u) + -j jtY A (t,u) 



which involves the four functions \lV (t , u), ^ A (t , u), X A (t , u) and Y A (t , u) as the 
independent three-particle distribution amplitudes. Due to this change, Eq. (11) must 
be also replaced by 



4>' + {t) - + (<A- J) (MO +2« 

= It C du{u- 1) (V A (t,u) + Y A (t,u)) . 
Jo 



' dMt) d4>-(ty 

dz 2 dz 2 



All other equations and the conclusions presented in the Letter remain unchanged. 
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